begin

  fi = addfile("../MCV/mcv_uv200hPa_global.nc", "r")
  u = fi->u
  v = fi->v

  lat = fi->lat
  lon = fi->lon
  dv = uv2dv_cfd(u, v, lat, lon, 1)
  vr = uv2vr_cfd(u, v, lat, lon, 1)

  copy_VarCoords(u, dv)
  copy_VarCoords(u, vr)
  ncdf = addfile("../divvor_global_mcv_ncl.nc", "c")
  ncdf->div = dv
  ncdf->vor = vr
  

end
